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coupled-cluster method 

S. G. Porsevi'2, M. S. Safronovai, and M. G. Kozlov^ 



(N 

o 

(N 



o\ 






O 

a.' 



> 



en 

o 



% 



^Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716, USA 
^Petersburg Nuclear Physics Institute, Gatchina, Leningrad District, 188300, Russia 

(Dated: February 27, 2013) 

The recently developed method combining the configuration interaction and the coupled-cluster 
method was demonstrated to provide accurate treatment of correlation corrections in divalent atomic 
systems [M. S. Safronova, M. G. Kozlov, and C. W. Clark, Phys. Rev. Lett. 107, 143006 (2011)]. 
We have extended this approach to the calculation of quadrupole polarizabilities 02 and applied 
it to evaluate 0:2 for the ground state of Mg and Mg-like Si^"*". Performing the calculations in 
three different approximations of increasing accuracy allowed us to place the upper bounds on the 
uncertainty of the final results. The recommended values 0:2 (3s^ ^^o) = 35.86(13) a.u. for Si^^ and 
Q!2(3s^ ^So) ~ 814(3) a.u. for Mg are estimated to be accurate to 0.37%. Differences in quadrupole 
polarizability contributions in neutral Mg and Si^^ ion are discussed. 

PACS numbers: 31.15.ap, 31.15.ac, 31. 15. am, 31.15.bw 



I. INTRODUCTION 

Atomic polarizabilities are important for a number of 
applications ranging from studies of fundamental sym- 
metries to development of ultra-precise atomic clocks as 
well as modeling of properties of chemical compounds. 
As a result, the study of polarizabilities has been of in- 
creasing importance in recent years. A number of both 
experimental and theoretical methods exist for determi- 
nation of El electric-dipole polarizabilities. We refer the 
reader to a recent review ^Ij] and references therein for the 
discussion of the applications of El polarizabilities and 
methodologies to determine these quantities in systems 
with a few valence electrons. 

The quadrupole polarizabilities have been a subject of 
far lesser number of studies, and few high-precision val- 
ues exist. The previous theoretical studies used the ex- 
ponentially correlated Gaussian functions 2], the pseu- 
dopotential methods Q-Qi coupled-cluster methods 0- 
0], and combined the configuration interaction and the 
many body perturbation theory (CI-hMBPT) fsj. Most 
of these studies have been limited to either very light or 
monovalent atomic systems. There are ever fewer exper- 
imental studies of the quadrupole polarizabilities owing 
to intrinsic difficulty in measuring this property. 

The quadrupole polarizability term arises in the effec- 
tive potential for polarization interactions between the 
core and the Rydberg electrons which allows determina- 
tion of 0:2 from the analysis of the fine-structure intervals 
of the high-L Rydberg states using a resonant excitation 
Stark ionization spectroscopy (RESIS) method 0. This 
approach has been used to determine the ground state 
quadrupole polarizabilities of Mg+ [10] . Ba+ tlUlSl' ^^'^ 
Th^+ [1^ [ij] . The lowest electric-quadrupole matrix el- 
ements involving the ground state, which generally give 
the dominant contribution to 0:2 have been determined 
using the RESIS method for the 6s — 5dj transitions in 
Ba+ [15] and 5/5/2 — 6dj transitions in Th-^^ [iq . 



The difficulty of determining the quadrupole polariz- 
ability from RESIS experiments that led to initial dis- 
agreement with theory in Ba+ was analyzed in [12]. The 
discrepancy between theory and experimental values was 
resolved in the same work [12]. The Mg+ experimen- 
tal a2 value from [10] was found to be in disagreement 
with both coupled-cluster and configuration interaction 
with semi-empirical core polarization potential (CICP) 
results , while both theoretical values are in excellent 
agreement with each other [5|. The revised polarization 
plot analysis of Ref. [13] data carried out in [g yielded 
much lower value, that is only slightly outside of the com- 
bined error estimates. The RESIS Th^+ a2 value is in 
agreement with theory predictions taking into account 
the uncertainties, but the theoretical calculations are of 
low accuracy for this ion. 

In summary, RESIS method can be used to determine 
the quadrupole polarizabilities to a good precision, but 
more benchmark cases of comparison with high-precision 
theory are needed. Therefore, it is important to develop 
approaches that are able to predict quadrupole polariz- 
abilities to high precision and evaluate the uncertainties 
of the final values. We note that quadrupole polarizabil- 
ity is a particularly good property for benchmark testing 
of theoretical approaches owing to generally large con- 
tributions of the correlation corrections. Accurate cal- 
culations and benchmark tests of divalent atomic sys- 
tem properties are of particular interest owing to appli- 
cations to atomic clock research 17], fundamental sym- 
metry studies 18], and quantum information 19]. 

The method combining the configuration interaction 
(CI) and the all-order coupled-cluster method (CI-|-all- 
order) developed in Refs. [20, [21[ was applied to the 
calculation of the electric-dipole static polarizabilities 
and the corresponding blackbody radiation shifts in di- 
valent B+, A1+, In+, and T1+ in J22, 23]. We apphcd 
the same approach to calculate electric-dipole polariz- 
abilities of the several low-lying states of Si^"*" ion in 



Ref. [2J| . These polarizabilities were calculated with un- 
precedented ~ 0.1% accuracies demonstrating a great po- 
tential and a high efficiency of the CI -I- all-order method 
in calculation of such quantities. Therefore, we select the 
ground state electric-quadrupole (E2) polarizabilities of 
Mg and Si^+ as the benchmark test cases for the approach 
developed in the present work, since the high precision of 
the theory is expected for these cases. We note that the 
method developed in this work can be used to calculate 
quadrupole polarizabilities of the low-lying states of any 
divalent or trivalent system. 

The E2 polarizability of the Mg ground state was calcu- 
lated using the CI-^MBPT method in Ref. @. The core- 
valence correlations were included explicitly in the second 
order of MBPT, while the higher orders of MBPT were 
included using the energy fitting described in Ref. [8|. 
The final resuh was 02 (Ss^ ^^o) = 812(6) a.u.. In this 
work, we carried out pure ab initio calculation of this 
quantity in the framework of the CI -I- all-order approach. 
The resulting ab initio value 814(3) a.u. is in an excel- 
lent agreement with the result reported earlier. A small 
difference between the central values can be attributed to 
the Breit interaction included in the present work, ab ini- 
tio treatment of the higher-order correlation corrections, 
the greater basis set, and the greater CI space. For Si^+, 
we obtained 02(85^ ^Sq) = 35.86(13) a.u. We will discuss 
the details of the method and evaluation of uncertainties 
in the following sections. 

To the best of our knowledge, there are no experi- 
mental data for the 02 (3s^ %) for Si^+ and Mg. At 
the same time the theoretical accuracy of calculation of 
these quantities achieved in this work is sufficiently high 
and, respectively, the E2 polarizabilities of the (3s^ ^Sq) 
state in Si^'^ and Mg present a good opportunity for a 
high-precision benchmark test of theory and experiment. 
Since Si^"*" electric-dipole polarizability was already accu- 
rately determined using the RESIS experiment [l2], it is 
a likely candidate for a benchmark test of determination 
of E2 polarizability by the RESIS method. Mg is also 
of particular interest since it is considered to be a good 
candidate for the development of the atomic clock. 



II. METHOD 

The details of the CH- all-order method were presented 
in [2ll . [23 . [24| . therefore, we provide only a brief outline 
of the approach. We start from the solution of the Dirac- 
Fock (DF) equations 

Hoipc = £ci>c, 

where Hq is the relativistic DF Hamiltonian |2l], ^^ and 
ipc and Ec are single-electron wave functions and energies. 
The dominant part of the Breit interaction was included 
in the self-consistency procedure [2^. 

The calculation was carried out in V^~^ approxima- 
tion, i.e. the self-consistent calculations were performed 
for the [ls^2s^2p^] closed core, and the 3s, 3p, 3d, 4s, 



4p, and M orbitals were formed in this potential. The 
B-spline basis set was formed in a spherical cavity with 
radius 60 a.u. and consisted of A^ = 35 orbitals for each 
partial wave up to Z = 5. The set of configurations was 
constructed by single and double excitations of the elec- 
trons from the main configuration 3s^ to the 4s — 23s, 
3p - 23p, 3d - 23d, 4/ - 23/, and 5g - 23g orbitals. Thus, 
five partial waves with the orbitals having the principal 
quantum numbers n < 23 were involved in the construc- 
tion of the CI space. We have verified that such CI space 
is numerically complete. 

The wave functions and the low-lying energy levels 
were determined by solving the multiparticle relativistic 
equation for two valence electrons [27| : 

i?cff(Sn)*n=-B„$„. (1) 

The effective Hamiltonian was defined as 

H,siE) = Hfc + ^{E), 

where iJpc is the Hamiltonian in the frozen-core approxi- 
mation. The energy- dependent operator S(-E) takes into 
account virtual core excitations. It is constructed us- 
ing the second order mariy-body perturbation theory in 
the CI-I-MBPT approach [23| and linearized coupled clus- 
ter single-double method in the CI-|-all-order approach 
|21| . The construction of the effective Hamiltonian in 
the CI-I-MBPT and CI-|-all-order approximations was de- 
scribed in detail in Refs. |2ll.l25|. 

The electric quadrupole polarizability 02 can be rep- 
resented in a general case as a sum of three parts 



a2 



vc 
"2 : 



(2) 



where a2 includes excitations of valence electrons, a2 is 
the ionic core polarizability, and a'!^'' is the small cor- 
rection to aj which subtracts out the excitations of core 
electrons into the occupied valence shells forbidden by the 
Pauli principle. In our case, a2'^ = because there are no 
nd orbitals in the core and, respectively, no excitations 
from the core to the 3s shell. The ionic core polariz- 
ability, Qfj, was evaluated in both the DF approximation 
and the RPA approximation. The difference between two 
these values can be used to estimate the uncertainty of 
this quantity. 

The static electric-quadrupole polarizability of the 
3s^ ^5*0 state can be written as 



a2 



2E 



(%|Qo|n)(n|Qo|%) 

En — En 



(3) 



In atomic units (m=?i=[e| = l) the electric quadrupole 



moment operator is determined as Qu 



C2.(n), 



where n = r/r and C2,y(n) are the normalized spheri- 
cal harmonics. 

The valence part of the polarizability, a2 , of the state 
|0) can be found by solving the inhomogeneous equation 
in the valence space, which is written as [28| 



(i?eff-So)|*)=Qcff|0) 



(4) 



and then calculating 



v^ = 2(0|(Qo)cff|*) 



(5) 



The effective quadrupole operator Qoff includes the 
random-phase approximation (RPA) corrections. 

Our analysis shows that the RPA corrections to values 
of q;2(3s^ ^5*0) for both Si^''" and Mg are very small (a few 
tenth of a percent). Therefore, all other smaller correc- 
tions to the effective operator including core-Brueckner, 
two-particle, structure radiation, and normalization cor- 
rections can be omitted at the present level of accuracy. 

We find that some caution is required in calculating the 
E2 polarizabilities by solving the inhomogeneous equa- 
tion in the valence space. First, the wave functions |^) 
and |0) are of the same parity. Second, the wave function 
|0) is the solution of the homogeneous equation 



(i/eff-So)|0)=0. 



(6) 



It is known that a general solution of an inhomoge- 
neous equation is a sum of a particular solution of the 
inhomogeneous equation and the general solution of the 
homogeneous equation. In our case the solution |^) of 
Eq. ([4]) is the sum of the particular solution of Eq. ^ 
which we denote as |^') and a solution of Eq. (jH]): 



\^) = \^')+f3\0), 



(7) 



where /3 is a numerical coefficient and |^') is assumed to 
be orthogonal to |0), i.e. (0|^') == 0. 

To separate out the particular solution |^') from the 
general solution, one needs to project |$) to the subspace 
orthogonal to |0) as 



I*') 



l*)-|0)(0|vl/). 



(8) 



We emphasize that such an admixture does not occur 
in calculating El polarizabilities, parity nonconserving 
amplitudes and other quantities for which JVP) and |0) 
are of opposite parity. For those operators, |\1/) and |0) 
belong to different subspaces from very beginning and 
automatically turn out to be orthogonal to each other 
(i.e.,/3 = 0). 

It seems that if |0) = \^So) we do not need to worry 
about the admixture of |0) to |^') because (0|((3o)cff|0) = 
and this admixture is removed from the final result. 
The problems arise because the factor /3 can be very large. 
In particular, in our case (3 was '^10^ which would lead 
to a numerical instability of the method and, finally, to 
a loss of accuracy in the straightforward implementation 
of the approach described above. 

We find a solution to this problem that can be imple- 
mented in the framework of our approach without addi- 
tional modifications of the method. The solution of the 
inhomogeneous equation, \'^), can be represented as a 
sum of projections to the states with definite total angu- 
lar momenta and written as [28] 



J+2 



w- E 



I* 



J',M), 



(9) 



J'=J,r 



where J and M are the total angular momentum and 
its projection of the state |0) and Jmin = max(0, J — 2). 
In our case, the only term in Eq. ^ that is of interest 
for us is |^j'=2,m)- All other terms do not contribute to 
Eq. ([5]). Thus, if we find the solutions of Eq. ([T]) belonging 
only to the subspace J' — 2,M — 2 we avoid the problem 
discusses above because the |*j'=2,m=2) and |0) = \^So) 
have different total angular momenta and cannot admix 
to each other. 



III. RESULTS 

First, we find the low- lying energy levels of Mg and 
Si^"^ . To estimate the accuracy of calculations we calcu- 
lated the energy levels in the CI, CI-I-MBPT, and Cl-t- all- 
order approximations. The results for Si^"*" were pre- 
sented in [2JI , where we demonstrated that the CI energy 
levels were already in good agreement with the experi- 
mental values. The maximum difference between the CI 
and experimental results did not exceed 2.2%. For Mg, 
the agreement with experiment at the CI level is only 
slightly worse. An inclusion of the core-valence corre- 
lations in the CI+MBPT and Cl+all-order calculations 
led to further substantial improvement of the theoretical 
energy levels. 

In Tables U and |TT1 we present the results obtained in 
the CI-I- all-order approximation for Si^+ and Mg, respec- 
tively, and compare them with experimental data. The 
two-electron binding energies are given in the first row of 
these tables, the energies in other rows are counted from 
the ground state. We find that the agreement between 
theoretical and experimental energy levels listed in these 
tables is extremely good, 0.05% or better for most of the 
levels. This is important for calculation of the quadrupole 
polarizability of the ground state because the low-lying 
levels give main contribution to this quantity. We note 
that the inclusion of the ng orbitals to the CI space is 
essential to obtain such high accuracy for all energy lev- 
els, including the singlet states, for Si^+. For Mg the ng 
orbitals can be omitted from the CI space with negligible 
loss of accuracy. 

In Table IIIII we list the contributions of several low- 
lying states to a2("^5'o) for Si^"*" and Mg. We also 
present the absolute values of the reduced matrix ele- 
ments |(^S'o||Q||n)|. For Si^"*", the two transitions from 
the ground state to the 3p^ ^D2 and 3s3d ^£'2 states con- 
tribute 97% of the final value. Such a large contribution 
of the 3p^ ^D2 state appears at first to be surprising be- 
cause the 3s^ ^Sq — 3p^ ^D2 is a two-electron transition. 
However, there is the large admixture (33% in probabil- 
ity) of the 3s3d configuration to the 3p^ configuration 
that explains such a large contribution of the 3p^ ^£'2 
state to a2(^<S'o). 

The breakdown of the E2 contributions is different for 
Mg, where the low-lying 3snd^D2 states with n = 3 — 5, 
listed in Table [1111 contribute 86% to a2(^S'o). For neu- 
tral Mg, whose electrons are more weakly bound to the 



TABLE I: Comparison between experimental [231 and theo- 
retical energy levels of Si^"^ in cm~^. Two-electron binding 
energies are given in the first row, energies in other rows are 
counted from the ground state. Results of the CI-|-all-order 
calculations are given in column labeled CI-l-All. Relative dif- 
ference of this calculation with experiment is given in the last 
column in %. 



State 
3s^ % 

Sp^'P2 

3s3d^D2 
3s3d^Di 
3s4s ^Si 
3p2 % 
3s4s % 
3s3d ^L»2 
3s3p^PS 
3s3p^Pi° 
SsSp^Pi 
3s3p iPi° 
3s4p^PS 
3s4p3Pi° 
3s4p^P2° 
3s4p ^Pi° 



Expt. 



CI-fAU 



634232 
122215 
129708 
129842 
130101 
142944 
142946 
142948 
153377 
153444 
159070 
165765 
52725 
52853 
53115 
82884 
175230 
175263 
175336 
176487 



634226 
122294 
129753 
129887 
130145 
142944 
142946 
142948 
153403 
153613 
159116 
165898 
52770 
52897 
53159 
82933 
175249 
175282 
175355 
176511 



Difference (%) 



-0.001% 
0.065% 
0.035% 
0.035% 
0.034% 
0.000% 
0.000% 
0.000% 
0.017% 
0.110% 
0.029% 
0.080% 
0.086% 
0.083% 
0.082% 
0.058% 
0.011% 
0.011% 
0.011% 
0.013% 



TABLE III: The contributions to the 3s^ % E2 polarizabili- 
ties (in a.u) in the CI -I- all-order approximation. The dom- 
inant contributions to the valence polarizabilities are listed 
separately with the corresponding absolute values of electric 
quadrupole reduced matrix elements given in columns labeled 
Q. The theoretical and experimental NIST transition ener- 
gies are given in columns A_Eth and A_Ecxp. The remaining 
contributions to 02(^^0) are given in rows Other. The con- 
tributions from the core are given in rows 02- The dominant 
contributions to 02 are calculated with the CI -I- all-order 
energies. 



Si' 



2+ 



Mg 



Contribution 



A£e: 



AEti 



a2 



3s' ^S-o - 


- 3p' 'D2 


122215 


122276 


5.200 


19.41 


3s' ^So - 


-3p"P2 


130137 


130101 


0.119 


0.01 


3s' % - 


-3s3d'-^D2 


143028 


142946 


0.002 


0.00 


3s' % - 


- 3s3d 'D2 


166247 


165765 


5.373 


15.24 


Other 










1.09 


a^2 










0.11 


Total 










35.86 


3s' % - 


-3s3d'D2 


46390 


46403 


18.65 


658.6 


3s' 'So - 


-3s3d^D2 


47939 


47957 


0.01 


0.0 


3s' 'So - 


- 3s4d 'D2 


53105 


53135 


5.05 


42.2 


3s' 'So - 


-3s5d'D2 


56289 


56308 


1.52 


3.6 


Other 










109.0 


ai 










0.5 


Total 










813.9 



TABLE II: Comparison between experimental [23] and the- 
oretical energy levels of Mg in cm~^. Two-electron binding 
energies are given in the first row, energies in other rows are 
counted from the ground state. Results of the CI-|-all-order 
calculations are given in column labeled CI+AU. Relative dif- 
ference of this calculation with experiment is given in the last 
column in %. 



State 


Expt. 


CIH-AU 


Difference (%) 


3s^ 'So 


182939 


188288 


0.03% 


3s4s ^Si 


41197 


41184 


0.03% 


3s4s 'So 


43503 


43491 


0.03% 


3s3d 'D2 


46403 


46388 


0.03% 


3s3d^D2 


47957 


47933 


0.05% 


SsSd'^Di 


47957 


47933 


0.05% 


3s3d^Di 


47957 


47933 


0.05% 


3s5s ^Si 


51873 


51854 


0.04% 


3s5s 'So 


52556 


52541 


0.03% 


3s4d 'D2 


52047 


53114 


0.04% 


3s3p^PS 


21850 


21849 


0.01% 


3s3p^Pi° 


21870 


21869 


0.01% 


3s3p^P2° 


21911 


21909 


0.01% 


3s3p 'P? 


35051 


35044 


0.02% 


3s4p^Po° 


47841 


47823 


0.04% 


3s4p3Pi° 


47844 


47826 


0.04% 


3s4p'P2° 


47851 


47833 


0.04% 


3s4p 'P? 


49347 


49328 


0.04% 



nucleus, the contribution of the high-lying discrete states 
as well as the continuum is larger than for doubly-ionized 
Si. Therefore, using a sum-over-states approach with a 
few low-lying contributions in the sum over intermediate 
states allows to obtain the ground state E2 polarizabil- 
ity for Mg only with an accuracy of about 15% percent. 
To obtain this quantity with a higher accuracy a more 
sophisticated approach like a solution of inhomogeneous 
equation used in this work is required. 

The E2 polarizabilities of the ground state obtained 
in different approximations are given for Si^"*" and Mg 
in Table IIVI This table illustrates that the core- valence 
correlations included in the second order of MBPT in 
the CI-I-MBPT approximation and in all orders in the 
CI-|-all-order approximation change the results obtained 
in the CI approximation by only a few per cent (by 4.5% 
for Si2+ and by 8% for Mg). 

Since we use the numerically complete basis set and 
the saturated CI space, we take into account the valence- 
valence interactions almost exactly. The main source of 
uncertainty arises from the core- valence correlations. We 
conservatively estimate this uncertainty as the difference 
between the CI-I-MBPT and CI-l-All results presented in 
Table Hvl The core contributions, 02, are very small for 
both Si2+ and Mg (0.3% for Si2+ and less than 0.1% 
for Mg) . Even if we estimate (very conservatively) their 
uncertainties at the level of 10%, their contribution to 
the uncertainty budget is negligible. Our final recom- 
mended values of q;2(3s^ ^Sq) are 35.86(13) a.u. for Si^"*" 



TABLE IV: Comparison of the present recommended values 
of the ground state static E2 polarizabilities (in a.u.) in Mg 
and Si^"*" with other calculations. First three rows give ab 
initio results for the valence polarizabihties Q2 calculated in 
the CI, CI+MBPT, and Cl+all-order approximations. In the 
row A(MBPT - All) the differences of the Cl+MBPT and 
Cl+all-order values are presented. The contributions from 
the core are given in the row a^. The "Total" are the values 
obtained as the sum of the CI+AU values and Q2 . 









Si^+ 




Mg 


CI 






37.58 




888.4 


CI+MBPT 






35.88 




817.5 


CI+AU 






35.75 




814.3 


A(MBPT - All) 






0.13 




3.2 


"2 






0.11 




0.52 


Total 






35.86 




813.9 


Recommended value 






35.86(13) 




814(3) 


Other works 






35.74(36)" 




812(6)' 
813.9'= 
809.3'' 

828'= 


■'Ref. ^; 'Ref. 


[8|; 


'=Ref. 


1; ''Ref. 3l|; 


"^Ref. 


M- 



ities and applied it to evaluate the static E2 polarizabil- 
ities of the ground 3s^ ^Sq state of the doubly-ionized Si 
and neutral Mg. Our recommended values are a2(^So) = 
35.86(13) a.u. for Si2+ and 814(3) a.u. for Mg. To the 
best of our knowledge, these are the most accurate values 
of these quantities obtained so far. They are in excellent 
agreement with the previous calculation of Porsev and 
Derevianko Q for Mg and with the theoretical results 
obtained by CICP method in Refs. 0,1131. 

The method developed in this work can be used to cal- 
culate quadrupole polarizabilities of the low-lying states 
of any divalent or trivalent system. We hope that 
the present work will stimulate experimental studies of 
quadrupole polarizabilities of divalent systems using the 
resonant excitation Stark ionization spectroscopy and 
other methods for benchmark test of theory and experi- 
ment. 



and 814(3) a.u. for Mg. Note that in both cases our re- 
sults are in excellent agreement with other most accurate 
results obtained by Mitroy in Ref. [30] for Si^'^ and by 
Mitroy and Bromley in Ref. [4] for Mg. 
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